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Upper bounds of the Hausdorff volume of scalar gradient field graphs are derived by means of 
geometric measure theory. The approach reproduces that scalar gradient fields along a mean imposed 
scalar gradient become space filling for sufficiently high values of Schmidt numbers Sc. The bounds 
are consistent with findings from recent high-resolution numerical experiments for 1 < Sc < 64, 
but too rough when compared with numerical simulations. A Reynolds number dependence of the 
bounds is found due to the additional scalar gradient stretching term in the equation of motion, 
^j- ". PACS: 47.27.Eq, 47.53.+n, 02.40.-k 
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Recent direct numerical simulations (DNS) suggested that the passive scalar mixing in a turbulent flow becomes more 
isotropic when the Schmidt number, Sc = v/k, is increased to values larger than unity, but the Taylor Reynolds 
number, R\, of the advecting turbulent flow is kept constant [1,2]. Here v is the kinematic viscosity of the fluid and k 

■ the diffusivity of the passive scalar field 9(x,t). The scalar was driven by a mean scalar gradient, G = Ge a , in both 
simulations which causes deviations from isotropy of the small-scale statistics. For increasing Sc, scalar filaments can 
be advected to ever finer scales which steepens up the local gradients, gi = d{9. Large gradients (or fronts) that are 
aligned with the mean, i.e. g a , occur already for Sc ~ 1 or even below and are associated with characteristic scalar 

■ structures, so-called ramps and cliffs (for further references, see [3]). A return to isotropic mixing is then thought as a 
growing compensation of those pronounced positive fronts by an increasing number of steep negative gradients with 
increasing Sc. Consequently, the probability density function of g a gets more symmetric tails and odd order derivative 
moments decay. 

Our understanding of this process for advection in Navier-Stokes flows is still incomplete and an investigation of 
geometric properties of scalar gradient fields was started therefore recently [4] . Fractal and multifractal properties of 
scalar gradient level sets and related quantities were studied there in a series of high-resolution DNS where, e.g., a 
higher degree of local isotropy was found to be related to a flater spectrum of generalized (multifractal) dimensions. 

On the theoretical side, concepts of geometric measure theory [5] were used sucessfully for scalars in turbulence 
[6,7] and relations to the scaling behavior of low-order structure functions were established. An extension of the 
framework discussed the dependence of the geometric properties of scalar level sets on Schmidt (or Prandtl) number 
t-H , and on spatial separation scales, respectively [8]. The purpose of the present brief report is to step in at this point 
and to extend the approach to scalar gradient fields. We will focus on the physics in the Batchelor regime of scalar 
turbulence where the advecting flow is in its viscous subrange [9]. Scales below the Kolmogorov dissipation scale of 
fluid turbulence, rj = (z^/e) 1 / 4 , but larger than the Batchelor scale tjb = rj/Sc 1 / 2 are considered, e is the energy 
dissipation rate of the flow. We will derive an upper bound for the scaling dimension of the Hausdorff volume of 
I^h . the scalar gradient field graph. This can give us an idea of how the gradients are distributed spatially and how this 
distribution depends on Sc. Our findings will be compared finally with high resolution numerical data for Schmidt 
, numbers between 1 and 64 to test the sensitivity of the derived scaling dimension bounds. 

The geometric measure theory generalizes concepts of differential geometry to non-smooth hypersurfaces embedded 
in an Euclidian space. The central object is the graph of the field under consideration, which is defined as 
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T r = {(x,g a )\x eB r d R 3 , g a = g a (x)} . (1) 

In case of a smooth field this is a three-dimensional hypersurface embedded in four-dimensional space. In Fig. 1 such 
graph is shown over a two-dimensional plane that is vertical to the direction of G. We observe that the hypersurface 
is strongly folded and rough for larger scales thus suggesting fractal properties. 

The Hausdorff dimension of such a graph, D h, (or more precisely the box counting dimension) is obtained from the 
scaling behavior of its Hausdorff volume H (T r ) with respect to scale r, the radius of balls B r of the covering [5,10], 

H{Y r )= [ Vl + r 2 |V5a| 2 d 3 x~r D " . (2) 
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The relative Hausdorff volume H(T r )/V r can be estimated as 
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|V5 a | 2 d 3 x, (3) 
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where the Cauchy-Schwartz inequality was used and V r = 4nr 3 /3. The scalar gradient field g a (x,t) is measured in 
units of its root mean square (rms) value, i.e. g a — g a /y (da)v- Here V = L 3 where L is the outer scale of turbulence. 
Similarily G = G/y/^j^. 

Progress is made now with a substitution of the gradient term in (3) by means of the advection-diffusion equation 
for the scalar gradient component, 

[d t + Uidi - nd 2 } g a + (d a Ui)gi = -d a u a G . (4) 
From the resulting second order balance one gets 

|V£fa| 2 = ^[(-«i5i + «^)^ 

- 2g a {d a ufg t + d a u a G)\ , (5) 

where summation is carried out over index i only, but not over a. The time derivative is already omitted, because we 
will discuss the statistically stationary case. With substitution (5), inequality (3) becomes 



H(T ' } <Ji + 



V r y Atttk J b 



+ ndf) ^ - g a (d a Ui)gi - g a d a u a G 



d 3 x . (6) 



We will consider now the four integrals under the square root separately and denote them by I\, I2, I3, and I4. From 
(6), with incompressibility, and by applying the Gauss theorem, it follows for I\ 

A r = 47rr 2 is the surface content of B r and the surface normal vector points toward the origin of B r . It is possible to 
add uo = u(xo), the velocity at the center of B r for which (uo)as r = 0. At this point, it has to be assumed that the 
fluctuations g 2 a are equally distributed over the sphere dB r in order to get § g 2 uq ■ dA = 0. The application of the 
Cauchy-Schwartz inequality results in 



h<^ A i Si^Jf ((u-uo)-n)^. (8) 
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The first square root is a scale resolved scalar gradient flatness which will be discussed later and denoted by F±{r) = 
(9a) B r / (9a)v ■ Note that the r-dependence comes in via coarse graining over balls of varying radius r. The second 
term stands for the second order longitudinal velocity structure function, which is S\\(r) = er 2 /(15v) in the viscous 
range, i.e. on scales around and below the Kolmogorov dissipation scale r\. We find 
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h < —\ r 1 . (9) 
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Integral I2 in (6) can be written as 

I 2 = ^J B [~ 9a V 2 g a + \Vg a \ 2 ]d 3 X . (10) 
Green's formula for scalars u(x) and i>(x), 

/ u(x)V 2 w(x)d 3 x = j> u(x)Vu(x)-dA 
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Vm(x) • Vw(x)d ci x, (11) 
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with w(x) = u(x) = g a (x) is taken. Substitution into (10) and application of the Cauchy-Schwartz inequality result in 

h = ^—l g Q |V.g a |n-dA 

47rr J dBr 



3r^F^)^(\W~g a \ 2 ) dBr , (12) 



2 



with F 2 (r) = (ga)dB r /{9a)v an d V.g Q = |V.g a | n. We assume that the volume average as well as the surface average 
give the same results over scales r due to homogeneity of turbulence. 

The third integral, 7*3, can be estimated by the Cauchy-Schwartz inequality to 
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/ g a (d a Ui)gid 3 x 



<-Y.^9l~9l) B A(d a u l Y) BT . (13) 



It is reasonable to assume that the fluctuations along the mean scalar gradient are the largest so that the mixed 
terms can be estimated by (g^gDsr _• {gV)B r — F^r) for i ^ a. The velocity gradient can be treated by the energy 
dissipation averaged over balls B r . Here we estimated Y^=i {{d a Ui) 2 ) B r < ( e )s r / I/ - The notation e r = (e)s r is used 
for the following and consequently 
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Similarily one can proceed for the integral I4 



H<-- S j^. (14) 



In summary we get 



H(T r ) 3r 2 eFAr) , — , Wl - ,„ ( r 2 FAr)e r r 2 G F 2 (r)e r 



(16) 



Unfortunately left with a couple of unknown terms in this expression. Further progress can be made only by the use of 
dimensional arguments. One can expect that the energy dissipation field fluctuates most strongly around scale 77 and 
large velocity gradients will be smoothed by the finite viscosity on scales below. Therefore it is reasonable to assume 

where (77/L) = (20/3) 3 / 4 i?~ 3/2 [11]. 7 is a (universal) scaling exponent varying between 71 and 72. Clearly, the 
minimum exponent, 71, will be the dominant one. 
The following term in I 2 , can be simplified to 

(\Vg a \ 2 ) dBr < -ka*) 2 W = ^fi , (18) 
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and F 2 (r) follows in lines with (17) to 

Fi(r) = * ^ < ^ * C-f) 5 ' 1 , (19) 

where ee >r = {eg)dB T is the coarse grained scalar dissipation rate. With 775/77 = Sc -1 / 2 one gets 

r 3 1 s - 1 
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F 2 {r) < 
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Again the scaling exponent 5 E [Si,5 2 ]. F 4 (r) can be estimated as follows 

(gt) Br {g 4 a )B rTPIT ^ (L^ 



F ^) = w 1 = j^f-W) < — F 4 (L) (21) 
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where F 4 (L) — (gt)v/(9a)v> (9a)B r ^ G 4 /i] B , and (g 4 )v — G 4 /L 4 . G is the magnitude of the mean scalar gradient. 
It is also assumed that scalar gradient moments are at maximum around the Batchelor scale. Thus follows 

27 

Mr) < 7^Sc 2 R\F A {L) . (22) 
v ; ~ 8000 v ' 



The term J 3 will be much larger than I4 if for the mean scalar gradient holds, G <C \J F4 (r) / F 2 (r) , as it will be the 
case for the comparison with numerical experiments. Finally, Ii will always be subdominant compared to 7 3 due to 
e r » e. 

All the estimates are inserted now into (16) and we assume a scaling relation for the relative Hausdorff volume, 
H(T r )/V r ~ r Dll ~ 3 . Additionally, all scales are expressed in units of Kolmogorov dissipation length, r — r/rj. The 
Hausdorff dimension of the graph over two-dimensional balls is given by the additive law, D' H = Dh — 1, [12] when 
assuming almost isotropic graphs, 

_ d in(g(r f yv f ) 

°H- 2+ dhif 



^ 2+ dln7 1 ^ 



1 + 3 



^ Sc^ f + i(71_1) ^ ^ f» . (23) 
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The expression probes the local slope of the scaling relation. The upper bound of D' H will be determined by the 
leading term under the square root at every scale f and we expect the result 2 < D' H < 3. 

The missing scalar derivative flatness factor F 4 (L) can be evaluated from [1] and present simulations. The latter 
are conducted in a homogeneously sheared flow in which the scalar field of constant gradient was allowed to evolve 
according to the advection-diffusion equation [4,13]. The mean scalar gradient is kept the same in all runs resulting in 
dimensionless G = 0.37 at Sc = 1 and 0.05 at Sc — 64 so that G <C y/ F± {r)/F 2 (r) is justified. As the inset of Fig. 2 
shows, the flatness factors for both simulations are ~ 10 for Sc > 1 and will be taken constant for the following. 

The two minimum scaling exponents, 71 and Si, can be taken from experiments where the whole multifractal 
spectrum for e$(x,t) and e(x,t) was measured, respectively. Meneveau and Sreenivasan [14] found for experimental 
data at different Reynolds numbers that 71 is about 1/4. Prasad et al. found a value of Si w 2/5 in a high-Schmidt 
number scalar mixing experiment at Sc « 2000 [15]. 

Figure 2 plots (23) for three different values of Sc. For scales accessible to the simulations, i.e, larger than tjb the 
upper bound is found to be always 3. While for the scalar case the crossover of D' H from 2 to 3 takes place around 
i]B [8], it is present here for by far smaller scales. It is just the fairly rough estimate for the additional scalar gradient 
stretching term (last term on the l.h.s. of (4)) that causes the crossover scale, f c , to be smaller than ijb- When inserting 
the numbers, one gets f c ~ S'c~ 6/ ' 5 i?^ 213/ ' 80 , at which the second term of (23) starts to dominate. To illustrate this 
effect, we kept 1\ only and plot the corresponding D' H as dashed lines in Fig. 2. The bounds are shifted then by an 
order of magnitude, but still rough due to estimate (22). 

For comparison, we conducted the relative Hausdorff volume of the scalar gradient field from numerical simulation 
data. Two-dimensional "balls" are then squares of a dyadic grid of sidelength 2~- 7 x 2ir. The corresponding local slope 
of the scaling dimension is added in Fig. 2 and a plateau with a non-integer dimension of D' H rts 2.3 can be observed 
for intermediate scales. This feature cannot be reproduced by the upper bounds because the f-dependence of all four 
/ terms has integer powers only. Even for passive scalars, non-integer Dh were obtained only for inertial range scaling 
of the velocity structure function of ~ r 2 / 3 (cf. [8]). 

To summarize, we have discussed the applicability of the geometric measure theory to scalar gradient fields for 
high-Sc mixing. The upper bounds on Dh are consistent with simulation results but rough. Additionally, they are 
insensitive to the particular choice of the threshold value of the level set and thus cannot capture slight differences 
between positive and negative level sets of same magnitude which were discussed in [4]. Interestingly, a dependence 
of the bounds on the Taylor-Reynolds number results. Batchelor's original model for the visocous-convective range 
was thought to be insensitive to the physics in the inertial range, i.e., on scales larger 77 [9]. A detailed investigation 
of that problem will be a part of future work and further improvement and extension of the bounds might be possible 
therefore. 

The numerical computations were carried out on the IBM Blue Horizon at the San Diego Supercomputer Center 
within the NPACI program of the US National Science Foundation which we wish to acknowledge. Comments by 
C. R. Doering, B. Eckhardt, and K. R. Sreenivasan are acknowledged. 
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FIG. 1. The graph of the scalar gradient field g a which is shown here over the x — z plane at fixed j/o and to for numerical 
data at Sc — 16. Grid resolution was N = 512. 
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FIG. 2. Scale dependent upper bound of the Hausdorff dimension of scalar gradient level sets for different Schmidt numbers. 
Value of R\ = 87 is taken from [4]. Solid lines are for Sc = 64 (thick line), 6.4, and 0.64 from left to right. Dashed lines are for 
bounds at the same Sc if integral I\ would dominate in (23) as in the scalar case [8]. The Batchelor scale for Sc = 64 is marked 
as a dotted vertical line. The panel contains also the DNS findings for the scale resolved D' H for Sc = 64 with a plateau at a 
value of about 2.3 . Inset: Flatness factor F4,{L) of the scalar derivative along the mean scalar gradient, g a , as a function of 
the Schmidt number. The dotted line stands for the mean of all data points and is drawn at 12. Present data are indicated by 
triangles and data from [1] by asterisks. 
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